Crystal-like high frequency phonons in the amorphous phases of solid water 
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The high frequency dynamics of low- (LDA) and high-density amorphous-ice (HDA) and of cubic 
ice (7 C ) has been measured by inelastic X-ray Scattering (IXS) in the 1-^15 nm~ momentum transfer 
(Q) range. Sharp phonon-like excitations are observed, and the longitudinal acoustic branch is 
identified up to Q = 8 nm -1 in LDA and I c and up to 5 nm" 1 in HDA. The narrow width of these 
excitations is in sharp contrast with the broad features observed in all amorphous systems studied 
so far. The "crystal-like" behavior of amorphous ices, therefore, implies a considerable reduction in 
the number of decay channels available to sound-like excitations which is assimilated to low local 
disorder. 

PACS numbers: 61.12.Ld, 63.50.+X, 64.70.Kb, 78.70.Ck 



Amorphous polymorphism, i.e. the existence of two or 
more amorphous states in the phase diagram of a chem- 
ical substance, has recently attracted wide interest from 
the scientific community JT|. In water, amorphous poly- 
morphism has received particular attention as it was as- 
sociated with a possible phase separation in the deeply 
undercooled liquid - this is suggested by the idea that 
the two amorphous water phases may be identified with 
the glassy forms of two liquid phases. In systems like wa- 
ter, the ergodic to non-ergodic transition from the under- 
cooled liquid to the glass cannot be studied continuously 
as a function of temperature due to homogeneous crystal- 
lization. The glassy nature of the amorphous ice phases 
can, therefore, only be established indirectly, and this has 
contributed to raise controversies on the exact nature of 
the phase-diagram of water, and on the inter-relations 
among the different stable and metastable phases of this 
molecule ||. 

The growing perception that the glassy state is a very 
general state of condensed matter justifies the many ex- 
perimental, theoretical and simulation studies performed 
on glasses. Among them, there is a growing interest in 
the study of the collective dynamical properties at high 
frequency, i.e. in the wavelength regime approaching the 
interparticle distances. Here one expects that either a 
continuum description, valid in the hydrodynamic limit, 
or the ordered-medium (phonon) picture developed for 
crystalline materials will fail as a consequence of the dis- 
order. The deviations can then be used to characterize 
and better understand the glassy state. More specifically, 
there is a number of dynamic signatures which are associ- 
ated with the glassy state Q] . At very low temperatures 
two-level states give rise to tunneling phenomena || . At 
higher temperatures disorder scattering and anharmonic- 
ities play an increasingly important role leading to excess 
densities of states (Boson peak) and relaxational phe- 
nomena In). Recently, the overdamping of sound waves 



at high frequency was proposed as a further criteria to 
characterize the glassy behavior 0. The study of the 
collective dynamics of amorphous ice at wavelengths ap- 
proaching the distance among water molecules may con- 
tribute, therefore, to the understanding of the differences 
between these two phases, and may help to shed light on 
such issues as their relation to stable crystalline phases 
or on the existence of two different liquid phases. 

In the present letter we present the excitation spec- 
trum of both amorphous forms of water, i.e. high-density 
(HDA) and low-density amorphous ice (LDA), as de- 
termined by high resolution inelastic X-ray scattering 
(IXS) at momentum transfers (Q) comparable to the in- 
verse of the inter-molecular distance. The experiment has 
been performed at the very high energy resolution ID16 
beamline at the European Synchrotron Radiation Facil- 
ity in Grenoble. The instrument energy resolution was 
set to 1.6 meV full width half maximum (FWHM) using 
21.748 keV incident photons. The detection system was 
made up of five independent analyzer systems, displaced 
among each other by a constant angular offset corre- 
sponding to a Q-spacing of 3 nm" 1 . The Q resolution was 
set to 0.4 nm -1 FWHM. The energy scans at constant Q- 
transfer took about 150', and each Q point was obtained 
by a typical averaging of 3-5 scans. The transverse di- 
mensions of the beam at the sample were 0.15 x 0.3 mm 2 . 
Further experimental details can be found elsewhere ||. 

The HDA sample was obtained by pressurizing hexag- 
onal D2O ice at 77 K || beyond lOkbar using 
the piston-cylinder apparatus described previously [flo|| . 
Keeping the sample always at liquid nitrogen tempera- 
ture, the metastable compound was retrieved from the 
cell at ambient pressure in mm-size chunks (p — 1.17 ± 
0.02 g/cm 3 ), and placed into a steel container with two 
diametrically opposed openings to allow for the passage 
of the incident and scattered X-ray beam. The effective 
powder sample thickness along the beam was « 15 mm, 
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matching well the X-ray photo-absorption length. Once 
filled, the container was transferred under a continuous 
helium flow from the liquid nitrogen bath onto the cold- 
finger of a precooled closed-cycle refrigerator. Above 
pa 90 K HDA transforms with a strongly temperature 
dependent rate to LDA (p = 0.94 ± 0.02 g/cm 3 ), which, 
in turns, converts itself into cubic ice I c at pa 140 K. HDA 
and LDA, after proper annealing, were both measured at 
60 K, while I c after annealing was measured at 80 K. The 
similarity of the temperatures allows a direct comparison 
of the scattering intensities. 

The purity of each phase was checked by measuring 
the respective static structure factors, S(Q), reported in 
Fig. 1 with a momentum resolution of 0.04 nm~ . Both 
HDA and LDA samples show the known static struc- 
ture factor |ll|] with no signs of Bragg-peaks and, there- 
fore, they are free of crystalline ice XII fH . There is a 
pronounced small-angle signal in HDA which disappears 
upon annealing to LDA. Within our energy resolution, 
its origin is purely elastic as it will be seen by the in- 
elastic scans reported in the following. As only moderate 
changes in the elastic small angle signal were observed 
with neutrons |13|, it is questionable whether the strong 
effect observed here is a real bulk property of the HDA 
sample or whether its origin stems from surface scatter- 
ing in the powder sample. 

A selection of inelastic spectra obtained in the range 
Q < 8nm _1 are reported in Fig. 2 and Fig. 3. In both 
amorphous ice phases and in the cubic crystalline phase 
there are well pronounced resonances. Increasing Q one 
initially observes a single feature with a marked disper- 
sion. At Q-values higher than 5 nm , the spectra be- 
come more complicated by the appearance of additional 
inelastic features. As X-rays, like neutrons, couple di- 
rectly to the longitudinal component of the density fluc- 
tuations, the dispersive excitation observed from the low- 
est Q-values can be readily identified with the longitudi- 
nal acoustic-like branches. In HDA these excitations can 
be clearly observed up to about 5 run -1 , while, in LDA 
and I c , they can be safely distinguished from the other 
feature(s) up to 10 nm -1 . 

In order to derive the energy position, KQ(Q), and the 
width, hT(Q), of the excitations, they have been fitted 
by a Lorenzian function convoluted with the experimen- 
tally determined instrumental resolution function. This 
resolution function corresponds to the central peak at 
E = in Fig. 2, that, as discussed previously, is due 
to the static disorder in the samples and to small angle 
scattering from the sample environment (see Fig. 2 and 
Ref. pi). The best fit to the inelastic X-ray scattering 
data shows that the excitation energy Ml(Q) scales lin- 
early with the momentum transfer Q in the limit of small 
momentum transfers. ^From the derived linear relation, 
Q(Q) — cQ, it is possible to obtain the longitudinal sound 
velocities chda = 3550±50m/s, clda = 3550±50m/s 
and a c = 3750 ± 50m/s, respectively. The differences 



among these values are small, and this is in agreement 
with the similar Debye level found in these three materi- 
als by neutron scattering experiments p3| . The determi- 
nation from the best fits of the resonance widths, KT(Q), 
is more involved, and is complicated by two facts: i) the 
absolute values for T(Q) are highly correlated with the 
background, and ii), particularly in HDA, the features 
appearing with increasing Q on the low-frequency side of 
the acoustic-like peaks interfere with the fit. Therefore, 
no systematic dependence of T(Q) on Q has been estab- 
lished. In both LDA and HDA one observes, however, 
that the resonances become broader with increasing Q- 
transfer. In any case care has to be taken in interpreting 
such broadening because the line shape of the resonances 
is unknown. For a single phonon it could be assumed 
Lorenzian reflecting the finite lifetime of the excitation. 
However, when taking a powder average the line shape 
acquires a non-analytic form even for the crystal due to 
the anisotropy within the multidimensional dispersion- 
sheet. 

Despite the difficulty to give full account of the ob- 
served broadening in the measured spectra, one can de- 
duce from Figs. 2 and 3 the following important obser- 
vation: As seen by the naked eye, the width of the reso- 
nances remains by far smaller than the excitation energy, 
hfl(Q), for all the considered Q-values. This result is very 
different from the observation made so far in other glassy 
systems. There, in fact, one observes an acoustic phonon- 
like resonance with a linear dispersion of Q(Q) vs Q and a 
quadratic dispersion of T(Q) vs Q up to a value Q = Q m . 
Q m is defined by the relation fl(Q m ) ~ r(Q m ). At Q 
larger than Q m it is no longer possible to observe well 
defined excitations and the inelastic part of the spectrum 
is, at the most, a broad and structureless background. 

Let us take the position, Qm, of the first sharp diffrac- 
tion peak as an indicator for the extent of structural 
correlations in the two amorphous phases. This allows 
us to define a pseudo Brillouin zone |14] which extends 
to Qm/2 ps 10 nm -1 in HDA, and to Q M /2 pa 8rrm _1 
in LDA. The quantity Qm/Qhi has so far always been 
found smaller than 0.5 [Q. In the two amorphous ice 
phases studied here, excitations are very well defined up 
to Q-transfer values approaching Qm, and, at least in 
LDA, the longitudinal acoustic-like branch can be iden- 
tified very well at least up to Q = 8nm _1 . In Ref. M, 
it has been suggested the possible existence of a relation 
between the value of Qm/Qu and the degree of fragility 
of the considered glass |{L5| . The high value of Q m /QM, 
which seems to approach unity in these two glasses, would 
imply that amorphous ices, and especially LDA, are the 
extreme end of fragile glasses. This is in contradiction 
to the view that water passes through an inflection in 
the deeply supercooled region where the liquid behavior 
changes from extremely fragile to strong M,[l7|. 

A further observation coming from the analysis of the 
spectra in Fig. 2 and 3 is that the spectra of LDA, at all 
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the considered Q-values, are very similar to those mea- 
sured in ice I c . In HDA the resemblance is less pro- 
nounced ]l9| . The similarity of LDA and I c holds not 
only for the resonances assigned to the acoustic mode, 
but also for the excitations at «7meV, which appear at 
higher Q and recall the transverse dynamics found in ice 
Ih and liquid water ||. In fact, as in these systems, the 
excitations set in around 5nm _1 and, as seen in Fig. 3, 
they become more pronounced beyond Qm- The trans- 
lational part of the density-of-states for Ih, I c and LDA, 
as obtained from INS spectra |f|, is peaked around 6.5 
meV (D2O). It is, therefore to be expected, that the INX 
spectra show excitations in this energy region at low Q 
due to Umklapp processes. These processes take place 
via the Bragg peaks in the crystalline state and via the 
static structure factor in the case of LDA [^0| . 

Despite the high definition of the low Q excitations in 
the two amorphous ice phases when compared to other 
glasses, one recovers a clear indication of the disordered 
character from the evolution of the inelastic spectrum at 
larger Q values. The spectra of LDA and HDA are less 
structured, lacking the sharp features observed in ice I c . 

To conclude, we reported on an IXS measurement of 
the S(Q, E) of the two known phases of amorphous ice. 
This has allowed to show that these two states of the wa- 
ter molecule possess a surprisingly crystal-like dynamic 
response. In both HDA and LDA the sound wave ex- 
citations are well defined. These experimental findings 
are in sharp contrast to the results found so far in other 
glasses, glass forming materials, liquids, dense gases and 
disordered materials in general. In these systems, in fact, 
an important broadening has always been observed in the 
inelastic part of the dynamic structure factor S(Q, E). In 
this Q-region the scattering experiment becomes sensi- 
tive to the topological disorder which opens decay chan- 
nels for sound excitations in addition to those available 
in the crystal. These channels are found practically ab- 
sent in the case of the amorphous ice phases indicating a 
very low degree of local disorder. Structural results are 
controversially discussed and to date do not give a clear 
picture of the topology J2l],|2^] ■ Experimentally only the 
pair correlation functions are directly accessible. Higher 
order correlations, among these the orientational corre- 
lation function, must be obtained in an indirect way |23j | 
e.g. via the dynamic response. We deduce from our data 
highly intact hydrogen bond networks both in LDA and 
to some lesser degree equally in HDA. Although an in- 
tact network in itself is no warranty for the absence of 
decay channels — e.g. an infinite random framework of 
corner-linked Si02-tetrahedra can undergo large phonon- 
assisted distortions — it seems a necessary condition. 
In LDA this network is perfectly annealed as it is not 
obtained via a fast quench from the liquid. Due to the 
constraints of the network the number of states the sys- 
tem can sample on the ps time scale should be small, i.e. 
there is a small configurational entropy, a view which is 



compatible with thermodynamic data p7|] . HDA is ex- 
pected to possess a larger configurational entropy, and 
on this basis one can justify that HDA has more "glassy 
behavior" than LDA. 

Apart from structure and bonding the dynamical prop- 
erties of the water molecule influence the decay of sound- 
like excitations both in crystalline and amorphous ice. 
We just want to point to the clear separation of trans- 
lational and librational bands which independent of the 
structural details arises from the very small moment of 
inertia of the water molecule. This separation closes de- 
cay channels — e.g. present in SiC>2 |2f| — involving res- 
onances of acoustic-like and librational modes, and may 
equally explain the absence of strong excess intensities 
in the inelastic neutron scattering data Ipjfl . Tempera- 
ture and the concomitant anharmonicities equally have 
to be given proper consideration in the discussion [p5| . 
In the end only detailed molecular dynamics calculations 
on well characterized ensembles combined with experi- 
ments on similar systems will be able to unambiguously 
give the reasons for the crystal like dynamic response of 
amorphous ice phases. 
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FIGURE CAPTIONS 



FIG. 1 - Static structure factor of high density (HDA) and low 
density (LDA) amorphous ice as measured on the inelastic 
X-ray beamline prior to the inelastic experiments. No signs 
of Bragg peaks are observed, indicating that both amor- 
phous phases are free of crystalline ice XII contaminations. 
Note the strong small-angle signal in HDA. The inset shows 
the diffraction pattern of cubic ice. 

FIG. 2 - Inelastic X-ray spectra of high density amorphous 
(HDA), low density amorphous (LDA) and crystalline cu- 
bic ice (I c ) at the indicated Q values lying mainly in the 
first pseudo Brillouin zone. The dashed line is a fit to the 
sinal using Lorcnzian lineshapes convoluted with the resolu- 
tion function (equally indicated). The solid line represents 
the inelastic contribution to the total fits. The numbers in 
brackets on the left of the elastic line give the elastic inten- 
sities in arbitrary units. Note the close resemblance of the 
sharp inelastic response of both amorphous phases to the 
crystalline phase. 

FIG. 3 - Inelastic X-ray spectra of high density amorphous 
(HDA), low density amorphous (LDA) and crystalline cu- 
bic ice (I c ) at high Q values in the second pseudo Brillouin 
zones. The dashed line gives the resolution function. The 
elastic intensities are given in brackets. For these high Q 
values the excitations become part of a broad intensity dis- 
tribution reminiscent of the density-of-states. 
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